Heavy baryons in the Skyrme model: the role of highly anharmonic collective motion 
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In the heavy quark and large N c limits, ordinary (exotic) heavy baryons can be considered as 
bound states of heavy mesons (anti-mesons)and chiral solitons. In these limits, the heavy mesons 
(or anti-mesons) and the chiral solitons are extremely heavy and are presumed to fall to the bottom 
of the effective potential. Previous studies have approximated the effective potential as harmonic 
about the minimum. However, when realistic masses for the heavy meson and chiral soliton are 
considered, longer range parts of the effective potential can become relevant. In this paper, we show 
that these longer-ranged effects yield effective wave functions which are qualitatively very different 
from those expected from the combined large iV c and heavy quark limits. These potentials can 
support bound heavy pentaquarks under some conditions. The consequence of these new energy 
states and wave functions on the Isgur-Wise function for the semileptonic decay of heavy baryons 
is also considered. 

PACS numbers: 12.39.Mk, 12.39.Hg, 12.39.Dc 



INTRODUCTION 



It has been known for a number of years that heavy 
baryons can be constructed in the heavy quark and large 
N c limits from the binding of heavy mesons with light 
chiral solitons 0, 0, El, 0, In these limits, both the 
heavy meson and the chiral soliton have large masses, 
and it is legitimate to describe the system in terms of 
a collective degree of freedom between the heavy meson 
and baryon. The large masses drive the particles to the 
bottom of the effective potential. In most cases consid- 
ered, this minimum occurs at the origin of the relative 
coordinate, causing both particles to be situated on top 
of one another. Previous attempts at describing heavy 
baryon physics within this model assumed that parti- 
cles would only experience small harmonic perturbations 
in the potential away from their minimum. The previ- 
ous works showed that the channels which constituted 
physical heavy baryons had attractive potentials while 
those channels which did not correspond to physical par- 
ticles had repulsive potentials 1, 2, 3]. Since the physical 
world does not correspond to the extreme heavy quark 
and large N c limits, which guided the previous work, it 
is interesting to ask the extent to which realistic masses 
drive the particles to the bottom of the potential (as ex- 
pected in the combined limit) when the harmonic ap- 
proximation to the effective potential is replaced by the 
leading order potential to all distances. In this paper, 
we look at a class of corrections to this picture. Similar 
issues have arisen for nucleon-nucleon forces in large N c 
[(|. We continue to assume that the heavy quark mass 
and N c are large enough to enable us to describe the dy- 
namics in terms of a collective degree of freedom which is 
describable as a non-relativistic effective potential. How- 
ever, we do not automatically assume that the masses 
are so large as to drive the particles into the vicinity of 
the the minimum of the potential. Instead, we focus on 
examining the consequences of considering the complete 



potential with realistic particle masses on heavy baryon 
physics. As we will show, the physics is qualitatively 
quite different from what one would expect if the world 
were close to the idealized limit. 

Baryons described as solitons in a chiral Lagrangian 
were first considered by Skyrme [7|. These chiral soli- 
tons have the correct quantum numbers as baryons in 
the large N c limit [8( when the Wess-Zumino term is in- 
cluded 0. While the extensive early calculations [1(3] 
focused on SU(2) solitons, the theory was also extended 
to include baryons with strange quarks. This was done 
by either considering the strange quark as light, and ex- 
tending the SU{2) chiral fields to SU{Z) [11 El, or by 
considering the strange quarks as heavy, compared to the 
lighter quarks, and strange baryons as bound states of the 
SU{2) soltion and K-mesons [lE[li|. 

To extend the theory to include even heavier quarks, 
such as charm and bottom, the former approach of ex- 
tending the soliton group structure is not reasonable since 
the mass of the heavy quarks is vastly different from the 
lighter three quarks. Therefore, the appropriate man- 
ner is the latter, put forth in the strange case by Callan 
and Klcbanov, where the heavy baryon states are consid- 
ered bound states of light quark chiral solitons and heavy 
mesons. For simplification, we will focus our attention on 
models made from SU(2) chiral solitons, as was done in 
previous work [3, 0, Furthermore, the theory must 
also exhibit heavy quark symmetry [TEl. HiL Il7l] . Thereby, 
the heavy meson is treated by heavy quark effective the- 
ory (HQET). HQET treats the heavy mesons as the dy- 
namical degrees of freedom and provides a systematic 
expansion, in powers of the hadronic scale, A, over the 
heavy quark mass, toq, to examine the Lagrangian and 
subsequent interactions. 

This method can also be used to bind heavy anti- 
mesons with chiral solitons to form pentaquark states. 
Such calculations have be en p erformed previously, for 
the case of strange quarks [l8[. This work showed that 
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the pentaquark channel was not sufficiently attractive 
to yield bound pentaquark states or prominent reso- 
nances. Previously, it has been argued on general model- 
independent grounds that in the extreme heavy quark 
and large N c limits, bound heavy pentaquark states must 
exist, primarily because the heavy particles fall to the 
bottom of the potential well However, it was shown 
that unlike in the extreme heavy quark and large N c 
limits in which the long-range one-pion-exchange poten- 
tial automatically binds heavy pentaquarks regardless of 
the details of the short distance interactions, for physical 
masses the existence of bound pentaquarks was highly 
sensitive to the details of the short distance interaction. 
The Skyrme type chiral soliton models provide one frame- 
work for treating the short distance physics, and thereby 
it is interesting to determine if such theories could sup- 
port heavy pentaquark binding. More significantly, the 
experience with the exotic channels suggests the possi- 
bility that the behavior in the non-exotic channels may 
also be substantially different from the limit on which the 
standard Skyrme analysis is based. 

The interest in determining the heavy baryon energy 
spectrum and wave function within a chiral soliton model 
has also led to studies of the Isgur-Wise function for tran- 
sitions between heavy baryons. Formally, the Isgur-Wise 
function describes the non-perturbative physics associ- 
ated with the form factor of the semilcptonic weak de- 
cay of Ab — > A c e~D e [2(|. Traditionally is is expressed 
as a function of the transition velocity. If the dynam- 
ics of the heavy baryons are known exactly, this function 
can be calculated for all velocities. The Isgur-Wise func- 
tion has been calculated previously when the effective 
potential was approximated as harmonic both in the con- 
text of the Skyrme model Q and in a model-independent 
context 0, HI . This paper considers the longer range ef- 
fects of the interaction between the heavy and light de- 
grees of freedom and leads to a non-universal form of 
the Isgur-Wise function dependent on the details of the 
interaction. 

We are by no means the first group to consider apply- 
ing the Skyrme model to attempt to explain either heavy 
baryons or exotic particles. In addition to 0,0, G|, heavy 
baryons spectroscopy within the Skyrme model has been 
previous considered by the work of Oh and his collabo- 
rators [21] . In recent years, a variety of work associated 
with pentaquarks states has been presented in the liter- 
ature. These works have included heavy pentaquarks, in 
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the context of several different models, 
and light pentaquark states either in the context of the 
Skyrme model, 
models, 



or in terms of other 
Exoctic dibaryon states have 
also been examined using the Skyrme model [HI- Ad- 
ditionally, decays of heavy baryons have been previously 
considered within the Skyrme model, 35]. 

The overall goal of this paper is to explore the proper- 
ties of heavy baryons in the Skyrme model in which the 
collective degree of freedom between the heavy meson 
and the remainder of the system is beyond the harmonic 
region of the ideal heavy quark and large N c limits. In 
the next section, we derive the complete effective poten- 
tial and show that when the full potential is considered, 
for realistic parameters, the heavy baryon wave function 
extends well beyond the region where the harmonic ap- 
proximation is applicable. This will be followed by a 
demonstration that the simplest type of interactions do 
not create bound states with realistic binding energies 
for the heavy baryons — indeed they give the same mass 
for ordinary heavy baryons and pentaquarks. This mo- 
tivates the inclusion of additional terms in the potential 
terms which allow for the correct binding energies and 
distinguish between heavy baryon and heavy pentaquark 
states. The consequences of these new potentials will be 
considered. Lastly, we will examine how these newly cal- 
culated wave functions for the heavy baryons influences 



the Isgur-Wise function. 



DERIVATION OF THE EFFECTIVE POTENTIAL 

The framework of this paper is based on the standard 
treatment of heavy baryons in the Skyrmion in terms of 
a collective degree of freedom between the heavy meson 
and the remainder of the system [H, 0, 0| . To begin the 
analysis of heavy baryons, the effective potential for this 
degree of freedom needs to be determined. The relevant 
Lagrangian for this purpose can be divided into the soli- 
ton sector and the HQET sector, 



L = L 



Skyrme 



c 



HQET- 



(1) 



The soliton portion determines the dynamics of the ordi- 
nary baryons. For concreteness we will consider the sim- 
plest such model; the one originally proposed by Skyrme: 



£ Skyrm c = ^ fl Tr^E^FJ] + _1_ Tr[(^S)Et, ( d „X)tf] 2 + ~ m \fl (Tr[E] - 2), 



(2) 



r 



where / w is the pion decay constant, e is the Skyrme pa- field Q- The last term provides the pion with a mass 
rameter, m v is the pion mass, and S is the chiral soliton 
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and fixes it to its physical value. As we will see below, 
this simplest model is not adequate phenomenologically. 
However, we will begin with an analysis of this model 
and consider more sophisticated models subsequently. 

Conventionally, the chiral soliton is treated classically 
as a first approximation and the ansatz taken for the form 
of the chiral field is, 



So (x) — Exp[ir • x F(r)], 



(3) 



where F(r) is the profile function determined by mini- 
mization of the soliton mass. Since this solution breaks 
both rotational and isorational invariance, any rotation 
within the isospin space will lead to a degenerate classical 
solution. That is, 



E(x) = A S (x) A~ 



(4) 



where A is an SU(2) matrix of the form A = ag + ir ■ a 
where ao and a (with a 2 , + a- a = 1) leads to a degenerate 
solution. This classical degeneracy is lifted when ao and 
a are promoted to quantum collective variables leading 
to baryons with the physical quantum numbers fioj]. 

However, there is a subtlety concerning the use of the 
quantization scheme of ANW for the regime considered 
here. There are two types of collective motion — the rota- 
tional/isorotational degrees of freedom of ANW and the 
collective vibration of the heavy meson off of the baryon. 
Provided that the natural time scales for these types of 
motion are very different they decouple. Formally, in the 
large N c and heavy quark limits, when the collective vi- 
brational degree of freedom resides near the bottom of 
the potential well and vibrates harmonically its natural 
time scale goes as A -1 / 2 [36| while the time scale of the 
ANW collective coordinate is A -1 , where A is an overall 
counting parameter for the combined limit: 



1 



M„ 
A 



A 



(5) 



where A is the characteristic hadron mass scale. For- 
mally, in the ideal limit, A — > and the two scales decou- 
ple. However, in the present study where we are explicitly 
looking for effects beyond the A — ► limit the situation 
is more problematic. 

As a practical matter, in order to proceed, we will as- 
sume that the motions do decouple. This can be justified 
in part on empirical grounds. The splitting between the 
S c and the A c in the extreme limit would be ascribed to 
a rotation excitation and hence of order A. The splitting 
between the first negative parity A c (the A c (2593)) and 
the A c is vibrational (in the ideal limit) and hence is of 



order A 1 / 2 which is perimetrically higher than the rota- 
tional excitation. In practice Me c — Ma c w 170 MeV 
while M Ac(2 593) - Ma c « 320 MeV. Thus, the ordering 
is what one expects. However, it is by no means obvious 
that one can legitimately take 320 MeV to be considered 
to be qualitatively large compared to 170 MeV. 

There is another reason why it is reasonable to treat 
the motions as though they decouple for certain qual- 
itative purposes. If our goal is to assess whether the 
standard Skyrmion treatment for charmed and bottom 
baryons (which implicitly works in the neighborhood of 
the combined limit) is self consistent, two issues arise. 
The first is whether the rotational and vibrational mo- 
tion decouple. The second is whether the collective vi- 
brational wave function is well localized near the bottom 
of the effective potential well. This paper is investigating 
the second question. As a logical matter, if one finds that 
the collective wave function does not remain well local- 
ized under the assumption that rotational and vibrational 
motion decouple, then there are important corrections to 
the standard treatment regardless of whether the assump- 
tion of decoupling is justified. Therefore, by showing that 
the collective wave function is not localized as in the ideal 
limits, our assumptions that the standard quantization 
method is applicable will be justified a posteriori. 

Having argued that the standard quantization of the 
collective coordinates is applicable for our purposes, we 
turn our attention to the heavy meson sector of the La- 
grangian. Near the heavy quark limit, the heavy mesons' 
momentum is essentially carried entirely by the heavy 
quark, and its properties are constrained by the emergent 
heavy quark symmetry [37J. Because of this, there is an 
approximate degeneracy between the spin-0 and spin-1 
heavy meson states. To incorporate this into the theory, 
HQET combines these two fields into one overall field, 



(i + jO 



(6) 



where w M is the four- velocity of the heavy quark with v 2 = 
1, P* M is the heavy vector meson field, P a is the heavy 
scalar meson field, and the index a is a light quark flavor 
index. The vector meson field is further constrained by 
v^P*^ — 0. Additionally, the field, H a , can be defined as 

H a = 7 °tft 7 ° = [Pftf + P a t 75 ]^. (7) 



Using these combined fields as the relevant degrees of 
freedom, the chiral Lagrangian for the heavy mesons can 
be written as, 
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£hqet = -iTrHaV^Ha + ^TiH a H b v» {pd^ba + ^TrH a H bl ^ 5 (£t0 M £) 



(8) 
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where the ellipsis denotes terms with more derivatives 
and inverse powers of the heavy quark mass, 1/mjj [HI, 
H, 17 1 . It should be noted that in this basis of the heavy 
meson, the heavy fields have an unusual transformation 
under parity Q. The field H a transforms under parity as 

H a {x°,x) -^^H^x^-xh^Ux^-x). (9) 

This transformation has the interesting property that 
when the soliton and the heavy meson are located at 
the same point, E^ = —1 in the transformation, while 
when they are separate E^ = 1. Thus, the heavy mesons 
act as though they have negative parity at long distances 
(as they must) but effectively as positive parity particles 
at short distances (when the anti-quark parity is also in- 
cluded) . 

The idealized A — > limit was uniform; i. e., the order- 
ing of the large N c and heavy quark limits was irrelevant 
|36|. It is hoped that in the current problem the ratio 
of the heavy quark mass to the ordinary baryon mass is 
equally irrelevant provided they are both large. If so, it 
is legitimate to calculate the effective potential assuming 
the heavy quark mass is infinite (but the nucleon mass is 
not). We will do that calculation here and subsequently 
verify that the result did not depend on this procedure. 
We will show the same effective potential arises if the 
nucleon mass is taken to be infinite and the heavy quark 
mass finite. 

Since we are considering the complete spatially ex- 
tended potential as a function of the relative separation 



distance between the particles, instead of the standard 
soliton ansatz, we will consider, 



E(£ - j?„(t)) = A(t)£ (2 - Mt))^ 1 (i), 
E (x - x (t)) = Exp[ir • — ^-r F(\x - x \)] ^ ^ 

\X — Xo\ 



where x is the position of the heavy meson, xq is the 
position of the soliton, and the soltion coordinate and 
the collective coordinates are time dependent while the 
heavy meson coordinate is not. This will fix the heavy 
meson to a specific location, which we will later choose 
to be the origin, and allows us to work in the rest frame 
of the heavy meson. Therefore the four-velocity of the 
heavy meson is v 1 * = (1,0). 

Even with this different choice for the soliton ansatz, 
the only term that produces an interaction between the 
heavy meson and the chiral soliton is 



(11) 



The summation is only over the spatial coordinates in the 
interaction as the spin trace with the temporal coordinate 
is zero. When the appropriate soliton ansatz is inserted 
into the interaction followed by some manipulation, the 
interaction term is written as 
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e jmk y*T 



m sin 2 (F) 



.4- 



ba 



(12) 



r 



where y — x — xq is the separation distance between the simplified by factoring out r from each term, 
soliton and the heavy meson. This can then be further 
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fl> = f <PZ r TrH a Hrf 1 ' i (AT i A- l ) ba 



y J y 



( 



F 



, sin(2F) \ sin(2F) 



) 



m 



■ h 3 + e Jlk y 



; *sin 2 (F) 



f /atoV (O. thaw (f (r- Sf2) + -^2 



k 3 + <'■'-.</■ 



sin 2 (F) 



(13) 
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By noting that the isospin operator of the heavy meson 
on the H field is, 



and the spin operator of the light degrees of freedom of 
the heavy meson on the H field is, 



P H 



■Hi 



(r j ) 



ba 



(14) 



S lH H b - -H b — , 



(15) 
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along with the fact that H^j 5 = in the rest 

frame of the H field, the interaction Hamiltonian can be 



These simplifications are similar to those performed in 
[H, whereas here we are considering the effective poten- 
tial to all distances. The integral can now be performed 
by explicitly fixing the heavy meson to be located at the 



The effective potential that we have just derived in 
Eq. (I17|) has three major aspects which are intercon- 
nected with each other. First, the potential is dependent 
on the spin and isospin of the light quarks in the heavy 
meson. Secondly, the term, Tr (At* A^t" 1 ), is related to 
the spin and isospin of the chiral soliton and is dependent 
on the collective coordinate quantization as well as the 
states being considered. The last part of the potential 
is the spatially dependent term. This term is a function 
of the separation distance, r, as well as the profile func- 
tion, F(r). The profile function can be derived numer- 
ically from the chiral soliton sector of the Lagrangian. 
Traditionally, it is achieved by minimizing the mass of 
the soliton either in the presence of a pion mass [38| or 
without a pion mass [loj . subject to the constraints that 
F(0) = —it and F(oo) = 0. Furthermore, if the pion 
wave function is expanded in powers of r and only terms 
of order r 2 are kept in the effective potential of Eq. (jTTj> , 
one can easily show that our potential reduces to the one 
considered by Jenkins, Manohar, and Wise, 0]. There- 
fore in the limit that the heavy meson and chiral soliton 
are close together, the effective potential is the same as 
previously considered. 



written as 



sin(2F) x sin(2F) sin%F) \ (w) 

2\y\ ) + 2\y\ V ffl )' { J 

I 

origin by equating TrH a H a and —S(x). Upon replacing 
the soliton position label, xq, with the more standard r, 
the effective potential reads: 



(17) 

I 

Let us turn our attention to the portion of the po- 
tential that is dependent on the chiral soliton, viz. 
Tr(^4rM _1 r m ). This term is dependent on the spin and 
isospin of the soliton, yet neither spin nor isospin are 
guaranteed to be valid quantum numbers of the opera- 
tor. That is, in some cases the chiral soliton term allows 
mixing between nucleon and delta states within the heavy 
baryon system. To simplify this issue, we will only con- 
sider the iso-scalar heavy baryons, as in the A#. This 
simplifies the problem because in order to create an iso- 
scalar from a heavy meson and a soliton, the soliton must 
have isospin- i. Additionally, large N c forces chiral soli- 
tons to exhibit the property that they have the same spin 
and isospin. Therefore, the only soliton that can bind 
with a heavy meson to form an iso-scalar heavy baryon 
has spin-i and isospin-^, or a nucleon. When the chi- 
ral soliton is confined to the nucleon sector, it can be 
shown that Tr^-rM-V" 1 ) is equivalent to -81% S^/3, 
where I™ and S % N are the isospin and spin of the nu- 
cleon, respectively. Making this replacement in the effec- 
tive potential and summing over repeated indices leads 
to a potential operator that reads: 



Hi = -gI%S{ H Tr(ArM _1 r m ) d 3 xTrH a H a 



y 3 y l 



ir 



I 



I 

V(r) = — —(Ih ■ In)[{Sih ■ r)(S N ■ r)(F ) + (Sih ■ S N ) — {Sm X S N ) ■ r . (18) 

6 V Zr Ir r J 



Having derived a potential operator, the problem re- operator. In order to determine the eigenstates of this 
duces to finding the eigenvalues and eigenstates of this 
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potential operator, let us consider states labelled by the the cross product term changes the spin state. Therefore, 



total isospin, /, the total spin, s, and the spin of the 
light degrees of freedom, s;. These states can be writ- 
ten as \I,s,Si). For a total isospin-0, we can construct 
three states; |0, ^,0), |0, \, 1), and |0, |, 1). From the po- 
tential, it is clear that total isospin is a good quantum 
number for the states, however, the spin of the light de- 
grees of freedom is not obviously a good choice here as 



instead of the simple state \I, s, s/), the appropriate wave 
function that should be considered has the form: 



*ao = [l + D{r)r- (S m x s" N )]\0, i 0)<f>(fl (19) 



for the |0, i, 0) state, and 



J 



tf ai = [1 + D'{r)r- (Sw x S~ N ) + E'{r){S m ■ f){S" N ■ f)]|0, -, l)</»(r) 



(20) 



r 



for the |0, \, 1) and |0, §, 1) states. There is a degeneracy 
for the light quark spin-1 states because of the degeneracy 
between the pseudo-scalar and vector heavy mesons in 
the heavy quark limit. It can be shown that the wave 
function for the light quark spin-0 state is in fact the 
eigenfunction of the potential operator when 



D(r)-. 

with an eigenvalue of 



-2(cos(F) + 1) 
r $m(F) 



2 r 



(21) 



(22) 



The effective potential for the isospin-0 light quark spin-1 
channel can be obtained in a similar manner. Here, the 
form given above is the eigenfunction when 

-4(1 + md 4 

rsin(i< ) r z 



with the eigenvalue 



rr / \ 9 jpu \ sin(F) 

Vai (r) = --F (r) - g . 

2 r 



(24) 



The previous discussion was based on taking the heavy 
meson mass to be arbitrarily large so that collective dy- 
namics involved the soliton moving. We have argued that 
the resulting dynamics ought to be independent of this 
assumption. To demonstrate this, the effective potential 
with the soliton's position held fixed can be calculated 
with the methods described above with a few caveats. 
First, since the heavy meson is now moving, the kinetic 
energy term of the heavy meson needs to be explicitly 
considered. From HQET [3^] the additional term in the 
Lagrangian for the heavy meson fields is 



^kinetic Tl_f/ a 



(Dl 



2M H 



(25) 



where is the covariant derivative perpendicular to 
the velocity and is defined as (-D^)ba = (D^)ba — v^{v ■ 



D)b a - The velocity of the heavy meson is given by v, 
and D M is the covariant derivative defined by (D^)b a = 
d^ha - |(S t 9' 1 S) 6o . The Roman indices are the light 
quark flavor indices, as before. Second, even though we 
are allowing the heavy meson to move, we would still 
want to be close to the heavy quark limit, therefore, the 
heavy meson's velocity will be small; w M = (1, e). 

The effective potential can still be derived from the in- 
teraction term as written in Eq. (jlip except the integral 
is now over xq — the soliton's position rather than the 
heavy meson's position. The use of only the spatial di- 
rections in this equation is still justified since corrections 
to this are of order e, which will remain small. The cal- 
culation proceeds as previously illustrated until Eq. (Ti"6|) . 
The substitution of the spin of the light quark in the 
heavy meson from the previous formula is still possible 
with again corrections of O(e). At this point, the integral 
can be performed analogously as before, but by fixing the 
soliton's position to be xq — 0. However, unlike before, 
we are left with the term TrHH in the expression of 
the effective potential. This term is dependent on the 
heavy meson wave function. However, the heavy meson 
wave function can be expressed as an exponential, i.e., 
H ~ Exp[i/(ir)], where f(x) is some unknown function of 
the heavy meson's position. With the appropriate nor- 
malization, we can thereby set TrHH = —1. Thus we 
have (practically) derived the same effective potential as 
in Eq. (TTT)) . The careful observer will notice that the 
potential derived with the soliton held fixed is identical 
with Eq. (fT7|) except for the sign of the last term which 
we will show will lead to the same physical system. 

The eigenfunctions and eigenvalues of the new poten- 
tial operator can be constructed just as before. How- 
ever, when we use the effective potential having held the 
soliton fixed instead, D{r) and D'{r) in Eqs. ([21) . ([23]) 
have the opposite sign. This sign difference compensates 
exactly the sign difference between the potentials dis- 
cussed above. Thereby both methods lead to the same 
physical effective potentials in Eqs. (|2"2")l . (f2"4"|) . Since both 
methods lead to the same physical effective potentials 
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we have demonstrated the commutativity of the large N c 
and heavy quark limits in this problem. 

To complete the discussion of these wave functions, we 
need to establish the correspondence between the \I,s,si) 
states and the physical Ah states. From the states' quan- 
tum numbers it is clear that the light quark spin-0 state, 
|0, |,0), corresponds to the ground state of Ah, while 
the light quark spin-1 states are spin excitations of this 
ground state which have yet to be observed. The ob- 
served excited states to A c , A+(2593) and A+(2625), 
would constitute radial excitations of the light quark 
spin-0 state. However, even with these apparently clear 
assignments, the parity of the heavy baryon states is not 
immediately clear in this language. The wave functions 
for both the light quark spin-0 and spin-1 states written 
above do not appear to have definite parity. Thus when 
the orbital momentum between the heavy meson and the 
soliton is considered in the I = state, the ground state 
wave function contains parts which are characteristic of 
both s- and p-wave states. The states achieve definite 
parity when we recall that the heavy meson state itself 
is negative under parity when near the soliton and pos- 
itive when far apart, as pointed out in Eq. ([9]) and the 
sentences following that equation. The wave functions 
in Eq. (fl"9|) and Eq. (|20|) show that when the heavy me- 
son and soliton are close together, the state has a pos- 
itive parity (positive from the s-wave, positive from the 
heavy meson) while when they are far apart, the state 
still has positive parity (negative from the p-wave, nega- 
tive from the heavy meson). Thereby, the ground states, 
and thus subsequent excitation, have the same parity as 
their physical particle states, which completes the assign- 
ment between wave functions and physical states. Note 
again, the subtlety associated with orbital momentum 
states since the orbital momentum is not a good quan- 
tum number. Henceforth in this paper we will label these 
states by the I used in the Schrodinger equation. This can 
be thought of as the orbital momentum state when the 
heavy meson and soliton are close together. In previous 
studies, since they were only concerned with small motion 
of the potential away from the heavy meson and soliton 
sitting at the same place, these long distance effects on 
the parity of the states did not matter. Thus, in their 
work the states could be clearly labelled by the orbital 
momentum between the heavy meson and the soliton. 

We have derived the effective potential of a heavy me- 
son and a nucleon in terms of the profile function, F(r), 
for the isospin-0 light quark spin-0 and spin-1 channels. 
Note that these potentials are completely radial. These 
potentials include short and long distance behavior for 
the binding which inherently has not been considered 
before. It should be noted that when both of these po- 
tentials are examined at short distances, they reduce to 
the potentials and values at the origin that have been 
previously identified [H, [3, Q . 
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FIG. 1: Calculated Skyrmion profile function, F(r). 

DETERMINATION OF BOUND STATES 

At this point, the effective potentials that we have con- 
structed can be used in a Schrodinger equation, and the 
bound states can be calculated. At the time of the pre- 
vious studies the heavy meson-soliton coupling, <?, was 
undetermined. In recent years this has been measured to 
be « 0.59 from the decays of D* meson into D mesons 
and pion emissions • We have assumed that this cou- 
pling is the same for B mesons (as it should be in the 
heavy quark limit) since an experimental determination 
via pion emission is not energetically possible. The phys- 
ical mass of the spin-0 heavy meson was used (1864 MeV 
for D meson and 5279 MeV for B meson [41(), while 
the mass of the soliton was calculated from the profile 
function. 

The short and long distance structure of the profile 
function, F(r), was obtained by examining the differen- 
tial equation that minimized the mass of the chiral soli- 
ton. From there, the profile function was constructed 
by parameterizing the functional form consistent with 
the short and long distance behavior with two param- 
eters. These two parameters where determined by an 
iterative method that minimized the mass of the soli- 
ton while keeping the nucleon-pion coupling, gA, and the 
pion decay constant, f^, constant. This procedure con- 
structed the Skyrmion profile function that is plotted in 
Fig. [Hand fixed the Skyrme parameter, e = 4.10, and the 
soliton mass, M = 949 MeV. 

The Schodinger equation with the appropriate effec- 
tive potential for this system was solved to observe the 
presence of bound states. Figure [2] shows the potential 
for the light quark spin-0 state with the harmonic os- 
cillator approximation overlaid. When the equation was 
solved, we found for the charm binding energy of 

155 MeV and for the bottom case, a binding energy of 
177 MeV. In both cases a weakly bound radial excited 
state was also observed; 6.18 MeV for charm and 19.32 
MeV for bottom. The observed ground states are more 
tightly bound than the ground state in the harmonic os- 
cillator approximation. Therefore the inclusion of the 
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entire potential increases the binding energy and favors 
a stronger bound state. Furthermore, the wave function 
of the nucleon in the ground state is much broader with 
the extended potential compared with the wave function 
of the harmonic oscillator (see Fig. [3]). This increase in 
the wave function breadth indicates that the nucleon is 
influenced by the long distance part of the potential. 

V [MeV] 

/ 

200 - / 
/ 

100 ■ / 




FIG. 2: Light quark spin-0 state effective potential (solid 
curve) with harmonic approximation (dashed curve) as a func- 
tion of separation distance. 

It is clear that for this particular model, both the shape 
of the wave function and the binding energy are vastly 
different when the entire effective potential is considered 
as compared to when the potential is assumed to be har- 
monic with the collective degree of freedom tightly lo- 
calized. Large amplitude motion clearly occurs and the 
standard analysis appears to be invalid for this system. 

For this system, orbital excited states can also be seen. 
Both the charm and bottom cases have an I = 1 excited 
state; the binding energy is 42.0 MeV (charm) and 68.0 
MeV (bottom). Neither system appears to have a bound 
I = 2 state. 

For the case when the light quark system carries spin 
1, no bound states were found. The potential in this 
channel has a strong repulsive core with a very shallow 
attractive region which appears to be too weak to support 
bound states. 

The previous calculations were performed using the as- 
sumption (valid in the heavy quark limit) that vector and 
pseudoscalar heavy mesons are degenerate. This is ob- 
viously not true in the physical case, that is, the D and 
the D* or the B and the B* have different masses. The 
calculation can be extended to include the physical mass 
splittings between the heavy meson states. When these 
splittings are included, there are still bound Ah states, 
however, the binding is weaker. The binding energy of 
the ground state is reduced by 88.7 MeV for charm and 
35.9 MeV for bottom. 

To reiterate, these results clearly illustrate that the ef- 
fective potential is not strong enough to localize the col- 
lective variable in the harmonic region for these models. 
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(a) Charm Wave Function 

*b[r] 



2 - \ 




(b) Bottom Wave Function 

FIG. 3: (a) Calculated normalized wave function for the 
ground state A c from the complete potential (solid curve) and 
the harmonic approximation (dashed curve), (b) Calculated 
normalized wave function for the ground state A& from the 
complete potential (solid curve) and the harmonic approxi- 
mation (dashed curve). 

The underlying assumption that the system is sufficiently 
close to the large N c and heavy quark limits to use the 
harmonic approximation is not justified. Clearly it is im- 
portant to see whether the breakdown of the combined 
heavy quark and large N c limit is generic for realistic nu- 
cleon and heavy quark masses. This is particularly true 
since the model considered has serious phenomenologi- 
cal flaws, as will be discussed below. The key question 
is whether the standard treatment works for "realistic" 
models. 



TOWARDS THE EFFECTIVE POTENTIAL IN 
REALISTIC MODELS 

The model considered in the previous section is unsat- 
isfactory in terms of phenomenology In the first place 
the mass of the heavy baryons is well off from the empir- 
ical ones. The relevant issue is not how large the frac- 
tional error is for the mass since a large fraction of the 
mass is simply from the heavy quark itself. The relevant 
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issue is the fractional error in the binding energy — i.e., 
the difference between the mass of a nucleon plus a heavy 
meson from the mass of the heavy baryon mass. In the 
model considered above M n + Mjj — M\ c = 155 MeV. 
In nature it is approximately 520 MeV. This is certainly 
not very good. 

Moreover, the model considered above has the feature 
that the interaction is identical under the exchange of a 
heavy meson to a heavy anti- meson. That is, our ex- 
tended potential will give the same bound states for or- 
dinary heavy baryons as for heavy pentaquark states. 
Clearly this is unphysical. While heavy pentaquark 
states are known to exist in the extreme heavy quark and 
large iV c limits [19], there exists no symmetry of QCD 
in that limit which implies a degeneracy between pen- 
taquarks and ordinary heavy baryons. Moreover, strong- 
interaction-stable heavy pentaquarks have not been de- 
tected despite intensive searches, suggesting that for re- 
alistic masses they do not exist. 

If we wish to make a more realistic model it is nec- 
essary to include additional interaction terms between 
the heavy meson and the soliton which split the ordinary 
heavy baryon from the pentaquark. Such an interaction 
should be strong enough to give heavy baryons with ap- 
proximately the correct mass. Heretofore, we have only 
included an interaction term between the pion fields and 
the heavy quark which was lowest order in the chiral ex- 
pansion. However, in the soliton there is no chiral power 
counting and thus no necessity to restrict the interaction 
to this term. The spirit of model building in Skyrme type 
models is to include a small number of terms to make 
the problem tractable. Although there is no systematic 
power counting, the hope is that one can get qualita- 
tively sensible results by choosing coefficients for these 
which compromise between the various observables. The 
simplest model realistic enough to get the binding of the 
heavy baryons correct while pushing up the mass of the 
heavy pentaquarks above threshold will require one new 
interaction term. 

The simplest term that we can consider, which distin- 
guishes between interactions between heavy mesons and 
anti-mesons, is coupling the light quark baryon current 
to the heavy quark vector current. Note that if a heavy 
anti-meson binds with a nucleon, the heavy quark cur- 
rent will switch sign compared with the heavy meson case 
while the baryon current (associated with the nucleon) 
will not. Therefore we add the term 

£baryo„ = g'TrB a H a ^B" (26) 

to the Lagrangian, with the baryon current, -B M , given 
by the standard form, 

B ^ = T^Tr[(St^S)(St5 Q E)(Et^E)], (27) 

Z47T 

where the notation £0123 = — e 0123 = 1 was used. By in- 
serting the functional form for the solitons into this term, 



and working in the rest frame of the heavy (anti)meson, 
the interacting potential can be derived as : 

VWa(0 = 9'B°(r) = ^ h ^l F '{r). (28) 

The coupling constant, g' , and its relative sign are un- 
known, but the sign will be chosen such that this po- 
tential is attractive for heavy meson-nucleon interac- 
tions while repulsive for heavy anti-meson-nucleon inter- 
actions. In the spirit of this class of model, we will tune 
g' in order to get a reasonable mass for the heavy baryon 
(either A& or A c ). Of course, this procedure is quite ad 
hoc, but Skyrme type models always require some ad hoc 
procedure. One useful check on whether the model ob- 
tained is sensible is whether the the coupling obtained 
has a natural size, i.e., whether it is of 0(1). 

When the Schodinger equation with the combined po- 
tential was used to calculate bound states, the coupling 
constants needed to have bound states with physical 
binding energies were determined to be -3.27 for A c and 
-3.34 for Aft. These are natural in size. Furthermore, the 
two coupling constants obtained via fitting for the two 
types of heavy baryons are quite close to one another; 
they differ by only 2% suggesting that the procedure is 
robust for physical values. For simplicity, the heavy me- 
son mass splittings were not included in this calculation. 

Because of the lack of the heavy meson mass split- 
ting, we would expect that the first excited state to be 
a superposition of the two physically observed excited 
states, A+(2593) and A+(2625). It is not unbelievable to 
surmise that this combined state would weight the two 
states by the relative spins of the states. This would 
lead to a mass of the combined state of 2614 MeV or an 
excitation energy of 328 MeV above the ground state. 
With the coupling constants considered here, the excita- 
tion energy from the combined potential was found to be 
324 MeV and 316 MeV for the charm and bottom cases, 
respectively. This is in surprisingly good agreement with 
the expected excitation from the physically seen states. 

With the inclusion of the baryon current term to the 
potential, the light quark spin-1 state becomes bound for 
the coupling constants considered here. The binding en- 
ergies are weak,; 125 MeV for charm and 173 MeV for 
bottom, compared with their spin-0 counterparts, 521.6 
MeV and 593.4 MeV for the A c and Aft, respectively. 
This binding is driven solely by the size of the baryon 
current coupling. The presence of a bound state in the 
spin-1 channel does not invalidate nor should it limit 
our method. It is not unreasonable to believe that the 
strength of the additional potentials that is needed to be 
included to drive the binding energy down to physical 
values might also bind light quark spin-1 states. There- 
fore the lack of any observation of these states could be 
associated with the difficulty in detection and not in the 
lack of the presence of these states. 

A harmonic approximation can also be made for the 
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combined potential with the provided coupling constants. 
The energy levels for the approximation are again differ- 
ent from that of the combined potential and the ground 
state is bound weaker. An examination of the shape of 
the wave function for the combined potential compared 
with the harmonic wave functions shows again that the 
wave functions for the combined potential are broader 
than the harmonic wave functions and extend beyond 
the harmonic region. We again find ourselves not driven 
to the effective potential minimum by the realistic masses 
suggesting that the physical states are not close to the 
extreme large N c and heavy quark limits. These ob- 
servations confirm the results from the previous section 
that the long distance potential greatly influences heavy 
baryon formation. 

Unlike with our original potential, the new combined 
potential is sensitive to the difference between regular 
heavy baryons and heavy pentaquarks. We can exam- 
ine whether a bound pentaquark state is possible within 
the combined potential by flipping the sign of the cou- 
pling constant of the baryon current term. When this is 
performed, a very weakly bound pentaquark state can be 
found with a binding energy of 33.4 MeV for a charm type 
pentaquark, and 42.9 MeV for a bottom type pentaquark. 
These are rather small binding energies compared with 
the binding energies of the regular heavy baryons, and are 
small compared to the deeply bound pentaquark states 
suggested by the large N c and heavy quark limits [lj|. 
Furthermore, the potential for the heavy pentaquark is 
repulsive at short distance creating the attractive region 
and localization of the wave function away from the ori- 
gin. Therefore, even in the ground state, the pentaquark 
formed here has a fixed separation between the heavy 
anti-meson and the nucleon. However, due to the strong 
repulsion at short distance this is a very delicate state. 
If we turn off the potential which is identical for the 
heavy baryon and heavy pentaquark states, and tune the 
baryon current coupling such that heavy baryons are still 
bound with the appropriate binding energy, we find that 
the heavy pentaquark effective potential is completely re- 
pulsive and thereby no heavy pentaquark state is bound. 
The coupling constants needed to achieve this condition 
were -5.52 for the charm case and -5.49 for the bottom 
case. These couplings are again neither unreasonably 
large nor small. This indicates that this model is capa- 
ble of binding a heavy pentaquark state, but with minor 
changes it is just as reasonable not to support a bound 
heavy pentaquark state. The relative ease for the state 
to become unbound furthers the point in 19] that heavy 
pentaquark binding with physically reasonable parame- 
ters is subject to the dynamical details of the model. 

We have justified the inclusion of an additional in- 
teraction term based upon the coupling of the heavy 
quark current to the baryon current. This interaction 
was derived for all separation distances. The coupling 
of the interaction was determined by having the ground 



state binding energy correspond to the physical value. 
With this potential in place, our simple model gives 
phenomenologically reasonable results concerning ordi- 
nary heavy baryons while simultaneously being reason- 
ably consistent and inconsistent with a bound heavy pen- 
taquark state. 

CALCULATION OF THE ISGUR-WISE 
FUNCTION 

The Isgur-Wise function is a universal function in the 
heavy quark limit that describes the semileptonic decay 
of all heavy hadrons in an SU(2A# ) multiplet (where Njj 
is the number of heavy flavors). In our case we can con- 
sider heavy hadrons of the A type [2(| ■ Thus the process 
in question is Af, — > A c e~D e Previous work has shown 
that the Isgur-Wise function is completely determined in 
the combined large N c and heavy quark limit Q • In this 
limit the effective potential is purely harmonic in nature, 
and the Isgur-Wise function would be dependent on only 
one parameter, the harmonic "spring constant" 0, Hf. 
Therefore it has been pointed out that if the excitation 
energy of heavy baryons were measured, the spring con- 
stant would be fixed and thus Isgur-Wise function would 
be completely determined up to higher order corrections. 

Unfortunately, we have shown thus far that for realistic 
parameters, the system does not remain in the harmonic 
region and the expansion appears to break down. Of 
course, from the wave functions we have already calcu- 
lated, the Isgur-Wise function can also be calculated, and 
again we would expect distinct results from the harmonic 
oscillator case. 

The most interesting aspect of the Isgur-Wise func- 
tion near the combined limit is related to the fact that 
this entire universal function is dependent on one mea- 
surable parameter. We wish to test the reliability of this 
result when one deviates from the ideal limit, viz., for the 
real world. Ideally one could test this by using the em- 
pirical value of the excitation energy of the first excited 
state plus the assumption that one is near the combined 
large N c and heavy quark limits to compute the Isgur- 
Wise function; this ought then to be compared with the 
experimental Isgur-Wise function for an empirical test. 
However, the Isgur-Wise form factors have not been mea- 
sured. 

In the absence of such data, we can still get some idea 
of how robust is the prediction by using "realistic" mod- 
els. Since these models indicate that that system is quite 
anharmonic, we do not expect a close agreement of the 
Isgur-Wise function with the harmonic approximation 
based on the energy of the first excited state. Instead 
of considering the entire Isgur-Wise function, we will fo- 
cus our attention on the curvature at zero recoil, p, of 
the Isgur-Wise function. We chose to do this as p is pro- 
portional to the mean radius squared in the heavy quark 
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limit, and thereby provides a single number with which In the case of harmonic wave functions this reduces to 
to compare and does so in a physically transparent way. 
The Isgur-Wise function in momentum space from [3[ 

is 

■qo = Jd 3 p<j)* c (p + m N v')(j) b (P)- (29) 
I 



VHo(z) = 



2V2 t f b /S " 3/S 



He 



Exp 



r 



where z = 



m N \v\ 



(V^ + V^) 1/2 ' 



(30) 



rriN is the mass of the nucleon, k is the harmonic cou- 
pling, v is the transfer velocity, and /i c and p,b are the 
reduced mass of the heavy meson-nucleon system for ei- 
ther charm or bottom, respectively. We can also express 
the Isgur-Wise function in terms of position space wave 
function as 



7]o = Id x ip* c (x)ipb{x) 



\ v \\ x \ 

TOjv|^||u| 



(31) 



These expressions for the Isgur-Wise function lead to the 
following leading order expressions for the curvature at 
zero recoil, 



Pho 



d 2 VHO 



a 2 "( z 



0) 



(32) 



for the harmonic case, and 



d 2 ij 

dz 2 



(z = 0) 



Pb 



/Me 



d x ifj*(x)ipb(x)\x\ 



(33) 

for the spatial wave function case. 

We will calculate p in three different manners. First, 
p can be calculated by approximating the potential by 
a harmonic oscillator to generate a harmonic coupling, 
K. Equation (j32j) can then be used for this value of k, 
along with physical values of the mass variables, to cal- 
culate p. Second, the wave functions of the ground state 
of the complete potential can be used directly in Eq. (f3"3"|) 
to calculate p. Of course, this is the proper way to cal- 
culate p and the Isgur-Wise function since the complete 
potential is being considered. We would expect a vast 
difference between the first two methods since we have 
been observing differences in the wave function between 
the harmonic wave functions and the complete poten- 
tial wave functions. Lastly, even though we know that 
the energy levels calculated by using the entire effective 
potential are not harmonic in nature, we could use the 
excitation energy of the first excited state to calculate 
k as though it were derived from a harmonic oscillator, 
and then use this value of k to calculate p using Eq. 



That is, instead of deriving the harmonic coupling from 
an approximation to the actual potential, we are con- 
structing a new harmonic potential, which is unrelated 
to the original harmonic potential and to the complete 
potential, except that the energy gap between the ground 
state and the first excited state of this new potential and 
the complete potential are identical. Of course, the third 
method has no underlying theoretical basis unless the 
system is harmonic (in which case it will agree with the 
first two). It is useful to consider it, however, because it 
simulates theoretically the empirical procedure outlined 
above. 

When the curvature of the Isgur-Wise function at zero 
recoil, p, was calculated in these manners, using phys- 
ical mass of the nucleon and the heavy baryons, the 
following results were obtained. The normalized wave 
functions are shown in Fig. [4J remembering that p is 
proportional to the expectation value of r 2 . Method 1 
yielded -53.8 x lCT 6 MeV~ 3/2 when g' was chosen to 
have the appropriate coupling constant for the bottom 
sector. Method 2 yields -121 x lCT 6 MeV" 3/2 . Note 
that by using the complete wave functions, the curva- 
ture is a factor of 2 larger than the case of using har- 
monic oscillator wave functions. Lastly, Method 3, of 
using the excitation energy from the full potential and 
assuming it came from a harmonic oscillator potential, 
yields — 111 x 10 _6 MeV~ 3 / 2 with the bottom case cou- 
pling constant. This result is different from the Method 
1 as one might have expected, but is quite similar to 
the Method 2 (which gives the "correct" result). This is 
quite surprising since Method 3 can only be justified via 
the harmonic approximation which appears to be badly 
violated (as seen from the result of Method 1). Given 
the fact that Method 3 simulates the natural way to test 
the harmonic approximation empirically, it is important 
to test whether the success of Method 3 in getting close 
to the correct result is a mere numerical accident for this 
model or whether it is a robust feature. As we will see it 
is quite robust. 

In order to test the extent to whether Method 3 generi- 
cally reproduces the correct result of Method 2, one needs 
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FIG. 4: (a) Calculated normalized wave functions for the 
ground state A c from the complete potential (solid curve), 
the harmonic approximation (smaller dashed curve), and the 
harmonic approximation from the excitation energy (larger 
dashed curve), (b) Calculated normalized wave functions for 
the ground state A& form the complete potential (solid curve), 
the harmonic approximation (smaller dashed curve), and the 
harmonic approximation from the excitation energy (larger 
dashed curve). 



to consider a wide variety of models and compare the re- 
sults. Since the effective potential only depends on the 
Skyrme profile function we have considered a variety of 
profile functions (which we constructed on an ad hoc basis 
entirely for the purpose of testing the validity of Method 
3). The curvature at zero recoil of the Isgur-Wise func- 
tion was calculated for these effective potentials using 
the three methods detailed above. The effective poten- 
tials that were considered are shown in Fig. [5] while the 
results of the calculations of p are presented in Table 
1. One can clearly see that with all effective potentials 
considered, the curvature at zero recoil calculated with 
harmonic wave functions is different from the calculations 
performed with the other two methods. However, for all 
the potentials considered, the latter two methods provide 
similar results. Although not depicted in Table 1, when 
an effective potential with a global minimum away from 
the origin is considered, Methods 2 and 3 give different re- 



sults. This can be attributed to the wave function calcu- 
lated with the complete effective potential being peaked 
away from the origin, where Method 3 assumes that the 
wave function is still peaked at the origin. 
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FIG. 5: Variety of effective potentials: (a) solid curve; (b) 
long dashed curve; (c) short dashed curve; (d) dotted curve; 
(e) dot-dashed curved. 



Potential 


Method 1 


Method 2 


Method 3 


a (Original) 


-53.8 


-121 


-111 


b 


-83.3 


-175 


-171 


c 


-109 


-216 


-215 


d 


-101 


-142 


-139 


e 


-132 


-225 


-231 



TABLE I: The curvature at zero recoil, p, of the Isgur-Wise 
function calculated in three different manners in units of 
10 _6 MeV _3,/2 . The first line is the original profile function we 
considered. All calculations were performed with the coupling 
constant, g' , chosen to bind At with the appropriate binding 
energy. 

The fact that Method 3 works so well, even when the 
system is quite anharmonic, appears to have some im- 
portant consequences. In the first place, it means the 
prediction of p from the excited state energy (using the 
harmonic approximation) may be expected to hold rea- 
sonably well and, hence, one has some real predictive 
power even if the system is rather anharmonic. The con- 
verse of this appears to be that the degree to which p is 
accurately predicted via Method 3 is a poor test of the 
degree to which the system is harmonic. 

In fact, the situation is a bit more subtle than this. 
Qualitatively it is clear what is happening: the anhar- 
monic nature of the potential lowers the excitation en- 
ergy compared to that of a harmonic potential with the 
same curvature at the minimum. Fitting this excitation 
with a harmonic oscillator means that the fitted oscillator 
will have a smaller curvature (i.e., spring constant) than 
the actual spring constant at the minimum. This has 
the effect of spreading out the wave function compared 
to the harmonic approximation based on the true curva- 
ture which in turn means a larger value of p; this acts 
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to simulate the true wave function which is also wider 
than the naive harmonic result. Thus, generically, the 
sign of the effect of anharmonicity on the excitation en- 
ergy and p helps explain the viability of Method 3. The 
degree to which the method works quantitatively may 
still seem remarkable, however. The quantitative suc- 
cess is at least partially understandable analytically. It 
is straightforward to calculate the leading order effect of 
the anharmonicity on both the excitation energy of the 
lowest excited state As it happens, the shift in the 
excitation energy exactly compensates the shift in p to 
leading order in the anharmonicity, and the according 
inaccuracies due to using Method 3 only appear at next- 
to-next-to leading order. Thus, the system can be rather 
anharmonic and Method 3 can remain reasonably accu- 
rate. It is nevertheless remarkable how well Method 3 
appears to work since the wave functions appear to be 
qualitatively quite different from the harmonic ones of 
Method 1. 

In conclusion, we have constructed the effective poten- 
tial for the binding of the heavy (anti-)meson and nucleon 
to form regular and exotic heavy baryons for variations 
of the Skyrme model. We have demonstrated that this 
effective potential gives excitation energies and collec- 
tive wave functions which are qualitatively different from 
those obtained with a harmonic oscillator approximation 
to the potential: the wave functions with realistic parti- 
cle masses are not concentrated near the potential mini- 
mum, as expected from the large N c and heavy quark 
limits. This indicates that the masses are not heavy 
enough to have heavy baryons exhibit the properties of 
these limits. Calculations with the full effective potential 
showed that heavy pentaquark states are possible in this 
class of model but the existence of bound pentaquarks 
depends sensitively on the details of the model studies. 
We also showed that despite strong anharmonicities the 
description of the Isgur-Wise function derived from the 
harmonic approximation works remarkably well provided 
that the effective harmonic coupling derived from the first 
excitation energy is used. 
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